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I.  Introduction 


For  many  physical  systems  it  is  important  to  follow  quantities  which  are  not  easily 
represented  in  the  common  orthogonal  coordinate  systems.  In  studies  of  plasmas  either  in  the 
earth’s  magnetosphere  or  in  solar  flares,  the  major  component  of  the  plasma  motion  is  along 
the  magnetic  field  line  which  is  approximately  dipolar.  These  are  examples  where  resolving  the 
flow  in  dipole  coordinates  has  two  major  benefits.  First,  by  using  the  most  natural  coordinate 
system,  it  preserves  the  physical  intuition  of  how  the  system  should  behave.  And  second,  use 
of  the  natural  coordinate  system  in  a numerical  simulation  will  minimize  the  numerical 
diffusion  due  to  interpolating  large  components  onto  other  coordinates. 


In  the  following  text  the  authors  have  compiled  the  results  of  their  derivations  for  the 
most  common  mathematical  formulas  and  operations  used  in  applications  of  a dipole  coordinate 
system.  The  results  are  given  for  a right-handed  coordinate  system  (V,L,<f>).  Previously  [1] 
the  coordinates  have  been  given  in  a slightly  different  order  (L,  V,  <f>),  where  this  latter  system 
is  left-handed.  We  note  that  the  results  for  all  of  the  operations  described  in  the  text  are  the 
same  in  either  system. 

II.  Geometry  of  the  Dipole  Coordinate  System 


Figure  1 shows  the  representation  of  a point  outside  a sphere  of  radius  r„  in  both  spherical 
( r , 0,  <f>)  and  dipole  ( V,  L,  <t>)  coordinates  where  the  relationship  between  the  coordinates 
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ile tines  a right-handed,  orthogonal,  curvlinear  dipole  system  with  coordinates  in  (he  range 
-1  < F < I,  1<  «.  0<  <t>  < 2tr 

By  relating  the  unit  vectors  of  the  two  systems  an  arbitary  vector  can  be  transformed  from  one 
system  to  the  other 
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where  8 - \fl  + 3cos20.  This  dipole  coordinate  system  has  a hybrid  representation  where  all 
coefficients  are  given  in  terms  of  spherical  quantities.  For  many  applications  this  is  the  sim- 
plest and  most  convenient  representation. 


A.  Metric  Coefficients  for  the  Dipole  System 
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B.  The  Differential  Arc  Length,  Area,  and  Volume  Elements. 
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C.  The  Derivatives  of  the  Coordinates  of  One  System 
with  Respect  to  the  Other  System 
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D.  Christoffel  Symbols 


Certain  vector  operations  are  simplified  by  the  introduction  of  these  symbols.  For  the 
dipole  coordinates,  there  are  four  independent  nonvanishing  components. 
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(1  + cos20). 
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III.  Vector  Operations  In  the  Dipole  Coordinate  System 

For  the  following  operations  / is  a scalar  function, 

/ - /( V.  L,  <t> ). 

A and  B are  vector  functions, 

(s)  - <-■  **> + {»;)^-  l-  *>!<-  * (»:!<*'  l *>i- 

and  T is  a corresponding  tensor  function, 

T-  T(V.  L.  <*>). 
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A.  Divergence  of  a Vector 
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C.  L.aplacian  of  a Scalar 
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I).  Curl  of  a Vector 


(V  x 4),  - ~±~-  JL  ( r sin#.4j ~±- 

/r.sm4tf  oL  /•  sin#  d</> 

(V  x A), -£-(r  sin#.4  J 

r sin#  9<2>  r4sin#  d I v 

(V  x 4)*  - - ' JL  siIL^  4 §L  l!  4 

r sin’#  91  ft  r„r’sin’#  BL  ft  ' 


K.  l.aplacian  of  a Vector 
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F.  Directional  Derivative  of  a Vector 
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G.  Divergence  of  a Tensor 
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IV.  Conclusion 


Most  of  these  expressions  have  been  used  and  tested  in  large  numerical  simulations  of  the 
ionosphere  and  magnetosphere.  We  hope  that  by  writing  the  operations  out  in  manual  form 
the  tedious  job  of  re-deriving  them  can  be  avoided  in  the  future. 
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